Kondo-temperature dependence of the Kondo splitting in a single-electron transistor 
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A Kondo peak in the differential conductance of a single-electron transistor is measured as a 
function of both magnetic field and the Kondo temperature. We observe that the Kondo splitting 
decreases logarithmically with Kondo temperature and that there exists a critical magnetic field 
Bq below which the Kondo peak does not split, in qualitative agreement with theory. However, we 
find that the magnitude of the prefactor of the logarithm is larger than predicted and is independent 
of B, in contradiction with theory. Our measurements also suggest that the value of Be is smaller 
than predicted. 

PACS numbers: 73.23.Hk, 72.15.Qm, 75.20.Hr 
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The many-electron Kondo state is formed when con- 
duction electrons screen a magnetic impurity. An im- 
portant tool for studying the Kondo effect is the single- 
electron transistor (SET)i*2ii4 which consists of a con- 
fined droplet of electrons, called an artificial atom or a 
quantum dot, coupled by tunnel barriers to two conduct- 
ing leads, called the source and the drain. When the 
quantum dot contains a net spin the conduction elec- 
trons in the leads screen the spin on the dot and en- 
hance the conductance of the SET. By using the gate 
electrodes that define the artificial atom to tune its pa- 
rameters, Goldhaber-Gordon et alA have tested the pre- 
dictions of renormalization group calculations and scaling 
theory and have shown that the equilibrium properties 
of the SET are quantitatively described by the Anderson 
Hamiltonian. 

SETs also provide the unique possibility of exploring 
non-equilibrium Kondo phenomena by applying a DC 
voltage Vds between the drain and the source. A sen- 
sitive tes1p§iL2i2iiS of theories of non-equilibrium Kondo 
physics is the measurement of the splitting in a magnetic 
field B of the spin- 1/2 Kondo peak in differential con- 
ductance as a function of Vd s - Calculations by Meir et 
al& predict that in a magnetic field the Kondo peaks oc- 
cur at eVds = ±A, where A = \g\firjB is the Zeeman 
energy for spin splitting and [Ib = 58 fieV/T is the Bohr 
magneton. We define Ax/e as half the separation in Vds 
between the two peaks, so Meir et al. predict Ak = A. 
Early measurements by Cronenwett et alA agree with 
this prediction. However calculations of spectral func- 
tions predict new features in the Kondo splitting. Costi£ 
predicts that the screening of the spin on the quantum 
dot by the conduction electrons should cause the Kondo 
peak to split only above a critical magnetic field Bq, 
which depends on the Kondo temperature Tk, the en- 
ergy scale that describes the strength of the screening. 
Moore and Wen 8 identify the field-induced splitting in 
their spectral function with the splitting in differential 
conductance and predict that the screening should cause 
Ak < A at all fields and that Ak decreases logarithmi- 
cally with increasing Tk- Recent observations by Kogan 
et a/^i have confirmed the presence of the critical mag- 
netic field Bq- However the latter measurements, as well 



as those of Zumbiihl et alJ&, show that Ak > A at high 
fields, in disagreement with theory. 

Here we extend the work of Kogan et (//.— by report- 
ing measurements of Ak as a function of the Kondo 
temperature Tk, as well as B. These measurements ex- 
ploit the property of SETs that Tk varies continuously 
with gate voltage, in contrast with conventional Kondo 
systems, for which Tk is determined by chemistry. We 
find Ak > A for high fields at the lowest Kondo tem- 
peratures, as observed previously. We also observe that 
Ak decreases as 1ii(Tk,o/Tk) with increasing Tk, where 
Tk,o is a constant. This logarithmic decrease in Ak with 
increasing Tk agrees qualitatively with theory. However, 
theory predicts that the decrease should be proportional 
to \g\fiBB ln(TK,o/TK), whereas we find the prefactor of 
the logarithm to be larger than predicted and indepen- 
dent of B. Our measurements also show that there is a 
critical magnetic field Be for splitting the Kondo peak, 
as predicted. However, we find that the predicted values 
of Bq as a function of Tk appear to be too large. 

The SET we study is fabricated from a heterostruc- 
ture consisting of an undoped GaAs buffer, followed by 
a 15 nm layer of Alo.3Gao.7As (5-doped twice with a total 
of 10 13 cm~ 2 of Si, and finally a 5 nm GaAs cap. The 
two-dimensional electron gas (2DEG) formed at the Al- 
GaAs/GaAs interface has an electron density of 8.1 x 10 11 
cm~ 2 and a mobility of 10 5 cm 2 /Vs at 4.2 K. Although 
the density of the 2DEG is high, magneto-transport mea- 
surements show that only a single subband is occupied. 
Electron-beam lithography is used to define the gate elec- 
trode pattern shown in the inset of Fig. [2fb) . Applying 
a negative voltage to these electrodes depletes the 2DEG 
underneath them and forms an artificial atom of about 50 
electrons isolated by two tunnel barriers from the remain- 
ing 2DEG regions, the source and drain leads. The elec- 
trochemical potential of the dot, as well as the coupling 
between the dot and the leads, can be tuned by changing 
the voltages on the electrodes. The voltage on the gate 
electrode g is denoted V g . The SET is measured in a 
75 [aW Oxford Instruments dilution refrigerator with an 
8 T magnet and a lowest electron temperature of about 
100 mK. To minimize orbital effects we align the 2DEG 
parallel to the magnetic field to within a few degrees. 
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We measure the differential conductance dl/dVds utiliz- 
ing standard lock- in techniques. 

To compare Ak to A we need an accurate measure- 
ment of the Zeeman splitting, and hence \g\, in our SET. 
Such measurements were reported for this SET in refer- 
ence 11 and the results are summarized here. We mea- 
sured the Zeeman splitting in this SET using the tradi- 
tional method of electron addition spectroscopy and the 
new, more precise method of inelastic spin-flip cotunncl- 
ing. These measurements showed that A was linear with 
field up to 7.8 T and the slope was given by \g\ — 0.16. 
We also used inelastic cotunneling to measure the Zee- 
man splitting in an identical device in a different dilu- 
tion refrigerator with a lower electron temperature and 
a larger magnet. In this setup, the device was aligned 
parallel to the magnetic field to better than 0.5 degree. 
Measurements of this sample showed that A was linear 
with field up to a magnetic field of 14 T and that \g\ 
was consistent with the value measured in the other de- 
vice. We attribute the small g-factor to the penetration 
of the electron wavefunction into the AlGaAs. We be- 
lieve the large electric fields produced by the relatively 
high electron density enhances this effect. Measurements 
of the Kondo splitting in a magnetic field in both devices 
showed that Ak increased linearly with field with a slope 
given by \g\, however in both devices Ak > A by about 
10 ^eV. 

We observe the Kondo effect by tuning the gate volt- 
ages so that the artificial atom is strongly coupled to the 
leads. Figure^a) shows an example of a plot of dl/dVds 
as a function of Vds and V g for this device. Enhance- 
ments in the conductance at Vds = in the middle of 
every second Coulomb blockade valley are clearly visi- 
ble. We infer that the valleys with Kondo peaks are 
those with an odd number of electrons. A sharp peak 
in dl/dVds at Vds = 0, for V g between the Coulomb 
charging peaks where resonant tunneling is prohibited 
by energy and charge conservation, is the hallmark of 
the Kondo effect. The gate lithography makes this device 
very small, resulting in typical excited state energies on 
the order 300 — 400 fieV. The large excited state energies 
can be seen explicitly in Fig. (Ha) by the location of the 
orbital inelastic cotunneling thresholds in the even val- 
leys. It has been demonstrated in both AlGaAs/GaAs 
and carbon nanotube SETs that a perpendicular mag- 
netic fieldi2ii&i£ii£ or gate voltageiSiiLiS can reduce the 
energy of the excited state and lead to higher spin Kondo 
states. The large excitation energy in our device makes 
this transition unlikely, especially since our magnetic field 
is parallel to the 2DEG, and we see no evidence for this 
transition in our data. For the remainder of this paper 
we focus on studying spin-1/2 Kondo physics. 

An example of a peak at B = T is shown in Fig.^c). 
Applying a magnetic field causes the Kondo peak to split 
into two peaks above and below Vd s = 0. This is shown 
in Fig. njb) ; while Fig-QIc) shows dI/dVd s vs Vds at var- 
ious magnetic fields. To study the splitting of a Kondo 
peak in dl/dVds as a function of Tk, we take advantage 
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FIG. 1: (a) Examples of Coulomb blockade diamonds when 
the artificial atom is strongly coupled to the leads. The Kondo 
effect shows up as sharp resonances near Vds = in the odd 
valleys. The inelastic cotunneling thresholds in the even val- 
leys indicate the energies of the excited states. These data 
are at different electrode voltages and during a different cool- 
down of the device from those in the rest of the paper, (b) 
Evolution of a Kondo peak as a function of magnetic field. 
Note that the splitting does not begin until B = 2.4 T. (c) 
dl/dVds vs Vds for B = T and for other values of B marked 
by the dash-dot lines in (b). For B < 2.4 T there is only one 
Kondo peak; above the threshold it splits. The traces have 
been offset by 0.08 e 2 /h for clarity. 



of the fact that Tk varies across a Coulomb valley. The 
dependence of Tk on the energy eg of an unpaired elec- 
tron in the artificial atom referenced to the Fermi level 
is predicted to b o 5 ' 19 



T K = 



2k f 



■exp[7re {e + U)/TU] 



(1) 



Here T is the full width at half-maximum of the Coulomb 
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charging peaks and U is the charging energy necessary to 
add an electron to the artificial atom. Goldhaber-Gordon 
et al£ have shown that their measurements of Tk are 
fit well by this equation. Expressing Tk in terms of its 
value in the middle of the Coulomb valley 



T K = T K ,oexp[7r(Ae) 2 /r[/] 



(2) 



In this equation, Tk,o is the Kondo temperature in the 
middle of the Coulomb valley, and Ae is the difference in 
the energy from its value in the middle of the valley. Ae 
is related to the gate voltage by Ae = a g e(V g — V gj0 ) 
where a g is the ratio of the gate capacitance to the total 
capacitance and V gj o is the value of V g in the middle 
of the valley. We can thus re-write the argument of the 
exponential in Equation (J2J as x(V g — V gi o) 2 where \ = 
ira^e 2 /TU . The Kondo temperature is a minimum in the 
middle of the Coulomb valley and increases exponentially 
close to the Coulomb blockade peaks. Thus, by varying 
V g , we can study the variation in Ak as a function of Tk- 
If the Kondo temperature gives a logarithmic correction 
to the splitting as predicted^ then we expect Ak to vary 
quadratically in gate voltage. 

To make quantitative comparisons with theory we need 
to determine Tk for a range of V g . To determine Tk at 
a fixed gate voltage we follow Goldhaber-Gordon et al~ 
and measure the differential conductance at Vds = (de- 
noted G) as a function of temperature. We then fit these 
data to an empirical form that gives a good approxima- 
tion to numerical renormalization group results 20 to ob- 
tain the Kondo temperature Tk- The empirical form is 
given by 



G(T) = Go[l + (2 1 /--i)(T/T K ) 2 



(3) 



where s = 0.22 for the spin- 1/2 Kondo effect. This fit is 
accurate only if we have data for temperatures as low as 
T ~ O.ITk- The lowest temperature we can achieve is 
about 100 mK, so we can reliably measure Kondo temper- 
atures larger than approximately 1 K. The Kondo tem- 
perature in the middle of the Coulomb blockade valley is 
significantly lower than this, so we cannot measure it di- 
rectly. However, we can measure Tk near the Coulomb 
blockade peaks where the Kondo temperature is much 
higher. Fitting the latter data for Tk as a function of 
V g to Equation J5J) we can determine Tk,o and X- 

Figure [3 illustrates this procedure. Figure Ufa) shows 
a Coulomb charging peak and half the Coulomb valley 
at two different temperatures. As the temperature is 
raised, the suppression of the Kondo effect removes the 
energy renormalization, causing the peak to move toward 
the bare resonance^Si We measure G as a function of 
temperature at several values of V g , one of which is shown 
in Fig. [2(b), and fit these data to the empirical form given 
by Equation to determine Tk- Finally, the data 
points in Fig. [2(c) shows Tk as a function of V g for 
the range in which it can be determined with confidence. 
Fitting these data to Equation J2J , shown by the dotted 
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FIG. 2: (a) dI/dV ds vs V g at V ds = (denoted G) for 
a Coulomb charging peak and half a Coulomb valley at two 
different temperatures. The solid line is a fit to Equations 
[5] and [3] and is discussed in the text, (b) G as a function 
of temperature for the value of V g marked by the dashed 
vertical line in (a). Fitting these data to Equation J3J gives 
Tk = 1.57 ± 0.05 K. The inset shows an electron micrograph 
of an SET similar to the one we have studied, (c) The dots 
show Tk determined as in (b) as a function of V g . The dotted 
line is the result of fitting these data to Equation The 
dashed line shows Tk as a function of V g using the results 
of the fit shown in (a). The solid black line shows Tk for the 
parameter values \ — 0.020 (mV) -2 and Tk,o = 250 mK. 



line in Fig. 12(c), we find T K ,o = 340 mK and x = 0.016 
(mVr 2 . 

Another way of determining Tk is to take advantage 
of the fact that Go is independent of gate voltage between 
the bare Coulomb resonances^. Combining Equations [2J 
and|3gives an equation for G(V g ) at a fixed temperature 
with fit parameters x,Tk,0j and Go- We use this equation 
to fit the G(V g ) data at T = 114 mK, shown in Fig. [2(a). 
The fit is shown as the solid line through the data and 
gives T K ,o = 145 mK and X = 0.026 (mV) -2 . A plot 
of Tk vs V g for these fit parameters is shown by the 
dashed line in Fig. [2(c). From these two methods, we 
arrive at the values x = 0.020±0.006 (mV)~ 2 and T K ,o = 
250 ±100 mK. A plot of T K vs V g determined from 
these parameters is shown as the solid line in Fig. [2(c). 
We can use our value of x to estimate T. From nearby 
Coulomb charging diamonds, we estimate U = 1.2 mV 
and ctg = 0.05. Combining this with our value of x — 
0.020 (mV) 2 we find that T « 330 /ieV. Our values are 
similar to those obtained by Goldhaber-Gordon et til.:' 
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FIG. 4: (a) Splitting of the Kondo peak in the center of the 
Coulomb valley as a function of magnetic field. The solid 
line shows A = \g\nsB for \g\ = 0.16 determined by inelastic 
cotunneling. Above ~4Twe find that Ak > \g\^,BB. Below 
B = 2.4 T no splitting is observed, (b) Curvature of Ak vs 
V g parabolas like that in Fig.|3Jb). The prediction of Moore 
and Wen 8 (dashed line) is that the curvature should decrease 
with field. 



FIG. 3: (a) A Kondo feature in a 7.5 T magnetic field. The 
Kondo peak at zero bias has clearly split into two peaks. One 
Coulomb charging peak is visible at V g = — 52 mV; the other 
peak is not visible, (b) Ak as a function of V g . The splitting 
is a maximum in the center of the valley, where the Kondo 
temperature is lowest. The gray line shows the result of fitting 
the data to a parabola. The dashed horizontal line indicates 
the Zeeman energy \q\hbB. Note that the splitting in the 
middle of the valley exceeds | g \fj, bB. (c) dl/dVds vs Vd B for 
the values of V g marked by the lines in (a). The traces 
have been shifted horizontally so that the peaks at negative 
Vds are aligned. The circles mark the positions of the peaks 
at positive Vds, showing that the splitting between the peaks 
changes with gate voltage. 



who found U = 1.9 ± 0.05 meV, a g = 0.069 ± 0.0015, 
and r = 280 ± 10 /xeV, and are consistent with having a 
somewhat larger artificial atom with stronger coupling to 
the leads. As expected for this stronger coupling, Tk,o m 
our SET is larger than that of Goldhaber-Gordon et al~ 
To study the dependence of Ak on Tk we measure 
dl/dVds as a function of Vds andV g at various magnetic 
fields. An example of such a measurement is shown in 
Fig-EIa), while Fig.^c) shows dI/dVd s vs Vds traces 
at the values of V g marked in Fig.^a). We determine 
the position of a peak by fitting the data points in the 
neighborhood of the peak to a parabola. The center of 
the parabola determines the peak position. At a given 



V g , Ak is half the separation between the positions of 
the positive and negative peaks. The values of Ak as a 
function of V g are shown in Fig. |3Jb) . These data are 
fit very well by a parabola, which is shown in the figure. 
As discussed above, a parabola is what one would expect 
from a correction to Ak that is logarithmic in the Kondo 
temperature. 

One important feature of these data is that for B > 4 
T, the maxima of the the Ak vs V g curves lie above 
A = \g\nsB, where the value of A is determined pre- 
cisely using inelastic cotunnelingiii By fitting Ak(V s ) 
to a parabola, we extract the maximum value of Ak, 
which is plotted as a function of magnetic field in Fig- 
ure Ufa). The solid line shows the value of A from co- 
tunneling measurements. These data show that for large 
magnetic fields, Ak > A by about 10 /ieV as observed 
previouslyiii*i£ 

Another important feature of these data is the cur- 
vature of the parabolas, which is plotted versus B in 
Fig. 01b). Moore and Wen£ predict that for 1 < 
A/fc B T K < 100, Ak/A = &ln(A/fc B T K ) + constant 
where & is a constant of proportionality and is approxi- 
mately 0.06. From Equation @ , Tk varies exponentially 
with gate voltage so that at fixed magnetic field the pre- 
diction is Ak = A K ,o - b(\g\(J, B B)x(Vg - V g , ) 2 , where 
Ak.o gives the field splitting at the center of the valley. 
This is the equation of a parabola in V g whose cur- 
vature becomes more negative with increasing magnetic 
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FIG. 5: (a) Drain-source sweeps in the middle of a Kondo 
valley at different magnetic fields. The onset of splitting is 
clear at B = 2.4 T but not visible at 2.2 T. The curves have 
been offset by O.ff e 2 /h for clarity, (b) Kondo valley with 
a 5 T field applied. The lines indicate the positions of the 
dl/dVds vs Vds traces shown in (c). (c) dl/dVds vs Vds at 
various gate voltages. In the center of the valley where Tk is 
a minimum the splitting is clear and symmetric. Closer to 
the charging peak the Kondo temperature is higher and the 
splitting is less pronounced, but still visible, (d) The dots 
mark the maximum Tk at which splitting is observed as a 
function of field. These data are extracted from data like 
those in (b) and (c) by measuring the gate voltage at which 
the splitting is barely visible. The conversion of this gate 
voltage to a Kondo temperature is explained in the text. The 
square data point at B — 2.4T marks the measurement of 
Be from the data in (a). Predictions for Be of Cost: 7 (solid 
line) and Moore and Wen 8 (dashed line) are shown. 



field. This prediction is plotted in Fig.^b) for x = 0.020 
(mV) 2 and clearly does not agree with the data. The 
magnitude of the curvature is much larger than predicted 
by theory, except at the highest fields, and does not grow 
with increasing magnetic field. The constant curvature of 
the data indicates that the correction to the Kondo split- 
ting is given by c1ii(Tk,o/Tk), where c is a constant; by 
comparing to the data we find c = 12 fieV. 

We can also use these data to test predictions about 
the threshold field Be necessary to split the Kondo 
peak. The existence of a threshold field is illustrated 
in Fig.^b). The peak does not start to split until B is 
greater than about 2.4 T. Figure[SJa) shows this thresh- 
old in more detail: at B = 2.0 T the peak has not yet 
split, but by 2.4 T, the splitting is clearly evident. We 
can also extract information about Be from data like 
those in Fig. |SJ For a given magnetic field there is some 



Kondo temperature Tk,c for which this field is the crit- 
ical field. For Tk below Tk,Cj the Kondo peak is split 
and for Tk greater than Tk,c the peak is not split. This 
is illustrated by the data in Fig. Efb) and in Fig. |3Jc), 
which shows some individual traces from Fig. |5fb). In 
the middle of the valley, where the Kondo temperature 
is lowest, the peak is clearly split. This splitting dis- 
appears toward the Coulomb charging peak where the 
Kondo temperature is higher. Because of the background 
from the Coulomb charging peak, it is difficult to reliably 
locate Tk,c- Instead, we identify a range of Kondo tem- 
peratures over which we are confident that the peak is 
split: this puts a lower bound on Tk,c- 

To extract this bound, we locate the most negative gate 
voltage at which we still confidently observe the Kondo 
splitting. For the data in Fig. E^b) and (c), this is at 
V g = —49.9 mV. To convert this to a Kondo temper- 
ature we use Equation (J2J. To find the center of the 
valley V gj o, we use the location of the maximum of the 
fit of Ak(V s ) to a parabola. We take into account the 
uncertainty in the determination of V gj rj to make the 
most conservative estimate of the lower bound. Knowing 
V gj o we can convert V g into a Kondo temperature using 
Equation J2J). The lower bounds on Tk.c found in this 
way using x = 0.020 (mV) -2 and Tk,o = 250 mK are 
shown as data points in Fig. OJd). The error bars show 
the Kondo temperatures extracted using the other two 
sets of x an d Tk,o shown in Fig. [j^c) and taking into 
account errors on these values. At Kondo temperatures 
below these data points, the Kondo peak is split. 

We can use our data to check predictions about Tk,c- 
Cost:' predicts that Be = O.SThwhm- where Thwhm is 
the half width at half-maximum of the Kondo resonance 
in the spectral function at T = and is relatedZ^ to 
Tk by Tk = 0.433Thwhm- This prediction is shown 
as the solid line in Fig. |3Jd): Costi predicts that be- 
low this line the Kondo peak should be split and above 
it the Kondo peak should not be split. However, this 
line lies below the lower bound set by some of our data 
points suggesting that at a given Tk, Be is smaller 
than predicted by theory. The dotted line in Fig. EJd) 
is based on the work of Moore and Wen, who predict 8 
the width and splitting of the spin-up component of the 
spectral function as a function of B/Tk- To model the 
Kondo peak arising from this spectral function we sum 
two lorentzians^ with the width and splitting given by 
the predictions of Moore and Wen. We find that the 
splitting appears at Be = Tk • This prediction is plotted 
as a dashed line in Fig. EJd) and lies below several data 
points at higher magnetic fields. 

The observation that Ak > A by Kogan et almL and 
Zumbiihl et alH demonstrates that we do not yet have 
a full theoretical understanding of the non-equilibrium 
Kondo effect. Our measurements of the Kondo tempera- 
ture dependence of Ak re-enforce this conclusion. While 
theory qualitatively describes the behavior of Ak with 
Tk, it does not provide an accurate quantitative de- 
scription. We find that Ak decreases logarithmically 
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with increasing Tk, but the prefactor of the logarithm is 
larger than expected at low fields and is independent of 
B. There is clear evidence of a critical field Bq for the 
splitting of the Kondo peak, again in qualitative agree- 
ment with theory. However, our data suggests that Be is 
smaller than predicted by theory. Our results provide ad- 
ditional challenges to the theory of the non-equilibrium 
Kondo effect. 
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